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Abstract 

We discuss the effects of noise on the timing and strength of superradiant Raman scattering from 
a small dense sample of atoms. We demonstrate a genuine quantum stochastic resonance effect, 
where the atomic response is largest for an appropriate quantum noise level. The peak scattering 
intensity per atom assumes its maximum for a specific non-zero value of quantum noise given by 
the square root of the number of atoms. 
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INTRODUCTION 



Stochastic resonance (SR) is a strongly surprising, yet very general effect in nonlinear 
dynamical systems. Against our naive understanding, the response of a system to an external 
driving can be facilitated if an appropriate amount of noise is added. In fact, the maximum 
of the response - the stochastic resonance - is found if the timescale of the noise matches 
an intrinsic time scale of the system By now, a stochastic resonance has been shown in 
a variety of classical systems ranging from neurosciences to geophysics. An overview can be 
found in the review articles 

In recent years there has been an increased interest in the emergence of stochastic res- 
onance phenomena in open quantum systems (see j^] for a review). These phenomena are 
quite diverse, as noise and dissipation can refer to different effects in quantum systems. Tem- 
perature induces a stochastic resonance in a bistable quantum system in a similar fashion as 
in the classical case: The response to a weak driving of the two states assumes its maximum 
for a finite value of the system's temperature jsl, Q]. Furthermore, temperature can also 
enhance genuine quantum features of an open system, as for example the entanglement in 
a quantum spin chain jsj. A different stochastic resonance effect can be found in terms of 
the coupling to the environment. In this case, the coherent response of a quantum system 
assumes its maximum for a non-zero coupling to the environment. This was demonstrated 

n 

for non-interacting quantum systems, e.g. in NMR experiments [9[, as well as for interacting 
many-body systems, as for example for the phase coherence and purity of a Bose-Einstein 



condensate 
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13|. 



In this article we will demonstrate another, even more subtle effect. We identify a stochas- 
tic resonance in the superradiant Raman scattering of a light pulse. Scattering is strongest 
if the time scale of the pulse is matched to the cooperative emission rate of the the sample, 
which is determined by the single-atom emission rate as well as the quantum shot noise. 
Thus we demonstrate a true quantum stochastic resonance, as the response of the sample 
assumes its maximum for a specific value of the quantum shot noise, that is, the uncertainty 
of the initial state. 
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THE DICKE MODEL 



Superradiance is a collective effect in a system of many atoms interacting via the light 
field. The spontaneous emission of a photon by one of the atoms stimulates the other atoms 
such that hght is emitted in a short pulse. The timescale of this pulse crucially depends on 
the presence of noise and dissipation, making superradiance experiments an ideal laboratory 
for stochastic effects in interacting many-body quantum systems. Furthermore, the basic 
principles of superradiance are generally valid and easy to understand, whereas the details 
of the dynamics are extremely involved [l^ . 

The essential features of superradiance are described by the seminal Dicke model intro- 



duced in 



15| (see 16| for a recent review). We consider light scattering from a small, dense 



sample of atoms, such that propagation effects of the light field can be neglected. Then 
one can eliminate the electric field and obtains a master equation for the emitters alone. 
Furthermore, the atoms are indistinguishable such that the dynamics can be described by 
collective variables. In particular, we define the collective atomic operators 

j 

j 

where \g)j and \e)j denote ground and excited state of the jth atom, respectively, and is 
the number of atoms in the sample. Then, J+ creates a symmetric excitation in the atomic 
sample and denotes the population imbalance between ground and excited state. The 
collective operators satisfy the commutation relations 

[J„4] = ±4, [J+,J-] = 2l. (2) 

For illustrative purposes it is often more suitable to work with the hermitian operators 

4 = ^(J+ + J_) and Jy = ^XJ+-J_). (3) 

The operators Jx,y,z form an angular momentum algebra with quantum number j = N/2. 
This representation is referred to as the Bloch representation. 

In terms of the collective atomic operators, the Dicke master equation for the dynamics 
of the atoms reads 



dt " 



P] - I {jJ-P + PJ+J- - 2Lph) . (4) 
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FIG. 1: Delay of the superradiant pulse as a function of the decay rate 7 (upper panels) and the 
atom number N (lower panels), (a) Intensity I{t) of the light pulse for = 20 atoms and three 
different values of 7. (b) Delay time ideiay as a function of 7 for = 20 atoms, (c) Intensity I{t) 
of the light pulse for 7 = 1 and three different values of N. (d) Delay time ideiay as a function of 
the atom number and 7 = 1. The dashed line is the mean- field prediction ([8]). 



and the intensity of the emitted hght is given by 



/(t)=7tr[J+J_p(t)]. 



(5) 



The decay rate 7 is a product of the spontaneous emission rate of a single atom and a 
geometrical factor. 

As the light intensity ([5]) is the essential observable measured in experiments, we briefly 
review some properties of the emitted superradiant pulse, in particular the time scale of 
the emission. Assuming that initially all atoms are in the excited state, each atom can 
spontaneously emit a photon. However, the emitted photons also stimulate the emission of 
further photons from the remaining excited atoms. Thus, the number of photons rapidly 
increases such that light is emitted in short intense pulse, as shown in Fig. [T] (a) and (c). 
The delay time tdeiay, at which the pulse assumes its maximum, obviously depends on the 
emission rate 7 as shown in Fig. [T] (b). 

Furthermore, the dynamics strongly depends on the quantum noise of the initial quantum 
state. When all atoms are in the excited state, the quantum uncertainty of the collective 



atomic variables are given by 

AJ^ = AJy = Vn. (6) 
The quantum shot noise, i.e. the relative uncertainty of the Bloch vector J, thus scales as 

lAJl 1 



(7) 



and vanishes in the classical limit — > oo. This noise drives the emission process such that 
the delay time of the superradiant pulse depends strongly on the atom number A^. As shown 
in Fig. [1] (d), the delay time first increases rapidly with A^ and then decreases again. The 
delay time can be inferred from a mean-field treatment of the atomic dynamics, taking into 
account the quantum fluctuations of the initial state ([6]). This approximation predicts that 



UeX.j = {lN)-HogN. (8) 

The mean- field prediction is compared to the numerically exact results in Fig. [1] (b) and (d), 
showing a good agreement when the atom number is large enough. 



DRIVEN SUPERRADIANT LIGHT SCATTERING 

The Dicke model shows a very rich dynamics when an external driving by a classical light 
field is included. The Hamiltonian describing the coherent driving is given by 

H = -Al + J-) , (9) 

where Q is the Rabi frequency and A is the detuning of the laser frequency with respect 
to the atomic transition frequency. If the driving is resonant (A = 0), the Dicke model 
undergoes a quantum phase transition in the thermodynamic limit N ^ oo. In this limit 
we can invoke a mean-field approximation, such that the dynamics of the expectation values 
Jz = {Jz) and J± = (J±) follows the equations of motion 

^ = -iAj+ + ir]j, + 7J+J, 

at 

^ = i^(J+-J-)-7^+^-, (10) 



5 




Rescaled emission rate 



FIG. 2: The steady state of the resonantly driven Dicke model, (a) Population imbalance {Jz)/N 
and (b) intensity of the emitted light ([5|) as a function of the rescaled emission rate 'jN/Q. Circles 
represent numerical results for = 20 atoms and = 1 , the dashed line is the mean- field prediction 
(fTTI) . The inset in (a) illustrates the system under consideration. 

On the resonance, the steady state is given by 





J. 



(11) 



for 2n < -fN 
N^/A for 2VL < -fN 

For strong driving or a weak decay, respectively, the population imbalance vanishes exactly. 
When the decay becomes stronger, > 2Q, a phase transition to an ordered state with 
(Jz) < takes place. This behavior is shown in Fig. |2] (a), where the population imbalance 
(Jz) is plotted as a function of the rescaled emission rate ^N/Vt. For finite particle number 
N , the transition becomes smeared out and (J^) is small but non-zero, even for 7A^ < 2Vt. 
The Dicke phase transition can be analyzed fully analytically for arbitrary particle numbers 



N 16|, so we do not go into detail here. An experimental observation of this phase transition 
has been reported only recently 17^ . 

In this article, we are primarily interested in the scattering light intensity, which can be 
viewed as the response of the atomic cloud to the coherent driving. Figure [2] (b) shows the 




FIG. 3: Superradiant Raman scattering by a sample of N A-type atoms, (a) Schematic represen- 
tation of the scattering process. The driving field 0,{t) is time-dependent with a pulse length T 
and maximum value 0,q. (b) The emitted superradiant pulse I{t) for three different values of 7 
and N = 20. (c) The emitted superradiant pulse I{t) for three different values of N and 7 = 1. 



intensity per atom as a function of the emission rate 7. One observes a clear maximum in 
the vicinity of the phase transition 7A^ ^ 2Q. This can be viewed as a stochastic resonance 
effect. The response of the quantum system assumes its maximum for a specific non-zero 
coupling to the environment, given by a specific value of 7. Moreover, the maximum is 
found if the timescale of the driving (2Q)~^ is matched to the timescale of the environment 
coupling (A^7)~^. 

We note that the scattered intensity essentially depends on the effective emission rate 
'jN/Q. Hence, the observed stochastic resonance is also present in the scattering of a single 
atom and when the atom number is varied. This effect is similar to the stochastic 
resonance observed in systems described by the Bloch equations, which are also governed by 

9HIII . Stochastic resonance effects in composite 



the competition of excitation and relaxation 



classical systems and their dependence on the number of constituents has been analyzed in 
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FIG. 4: Stochastic resonance in superradiant Raman scattering: The peak intensity per atom of 
the emitted pulse assmnes a maximum for a specific non-zero value of (a) the coupling rate to the 
environment 7 and (b) the magnitude of the quantum uncertainty ^/N . The remaining parameters 
are f^o = 1 (blue circles) and VLq = 2 (red squares), A = 1, T = 5 and (a) = 20 and (b) 7 = 1, 
respectively. Solid lines are shown to guide the eye. 
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STOCHASTIC RESONANCE IN PERIODICALLY DRIVEN SYSTEMS 

Superradiant light emission is essentially started by quantum fluctuations such that its 
timescale is determined by the quantum uncertainty A^^/^. Therefore we can identify another, 
much more subtle stochastic resonance effect. 

We consider superradiant Raman scattering of a light pulse by a dense sample of three- 
level atoms as shown in Fig. [3] (a). The atoms emit light collectively, when they decay 
from the excited state |e) to the ground state \g) as described by the Dicke master equation 
dl]). However, we assume that they are initially prepared in a metastable state which is 
coupled to the excited state by a classical light pulse, 

H=-^ (Jss + Jee) - ^ (Jse + ■ (12) 
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The collective operators used here are defined as usual 

Ja>, = ^J2\^)Ab\,, (13) 
i 

where a and b label one of the atomic states |s) or |e). For simplicity, we assume that the 
input pulse is given by fl(t) = f2o sin(7rt/T)^ in the interval t G [0,T]. We solve the master 
equation numerically using the Monte Carlo wave function method [20I averaging over 
500 (A^ < 15) resp. 100 (A^ > 20) random trajectories. 

Figure [3] (a) and (b) show the emitted light pulse for different values of 7 and the atom 
number N. In both cases, we observe that the peak intensity is largest for intermediate 
values of 7 or A^, respectively. Obviously, the peak intensity is rather small if the emission 
rate 7 and the atom number N are small. In this case the emission is too weak and most 
atoms remain in the initial state However, if 7 or are too large, then the emission 
starts 'too early', long before the driving Q{t) has reached its maximum. Also in this case, 
the peak intensity is limited. 

This is further analyzed in Fig. HI where the peak intensity per atom is plotted as a 
function of (a) the emission rate 7 and (6) the atom number A^. In both cases we observe 
a clear maximum when 7 or A^ assume intermediate values, such that the time scale of the 
superradiant emission matches the timescale of the driving. The second case is particularly 
interesting, as the maximum is found for a specific non-zero strength of the quantum un- 
certainty A^^/^. The observed stochastic resonance is thus a genuine quantum effect, as it is 
observed in terms of the quantum shot noise. 



CONCLUSION 



In this article we have analyzed different stochastic resonance effects in superradiant light 
scattering. We have identified a genuine quantum stochastic resonance, where the response 
of the atomic sample to an external driving assumes its maximum for a specific non-zero 
strength of the quantum shot noise N~^^'^. The response is quantified by the scattered light 
intensity per atom or its peak value, respectively. 

The mechanism generating the stochastic resonance is a matching of the timescales of 
the driving and the superradiant emission. The latter depends on the emission rate 7, but 
also on the quantum shot noise of the sample, as the the emission process is essentially 
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started by quantum fluctuations. Thus we can observe a maximum of the response, the 
stochastic resonance, in terms of both the couphng rate to the environment 7 and the 
quantum uncertainty 

Acknowledgement s 

I thank A Hilhard and J H Miiller for stimulating discussions. Financial support from the 
Deutsche Forschungsgemeinschaft (grant number DFG WI 3415/1-1) and the Max Planck 
society is gratefully acknowledged. 



[1] Benzi R, Sutera A and Vulpiani A (1981) J. Phys. A: Math. Gen. 14 L453 

[2] Wiesenfeld K and Moss F (1995) Nature (London) 373 33 

[3] Dykman M I, Luchinsky D G, Mannella R, McClintock P V E, Stein N D and Stocks N G 

(1995) Nuovo Cimento D 17 661 

[4] Gammaitoni L, Hanggi P, Jung P and Marchesoni F (1998) Rev. Mod. Phys. 70 223 

[5] Wellens T, Shatokhin V and Buchleitner A (2004) Rep. Prog. Phys. 67 45 

[6] Lofstedt R and Coppersmith S N (1994) Phys. Rev. Lett. 72 1947 

[7] Wellens T and Buchleitner A (2000) Phys. Rev. Lett. 84 5118 

[8] Huelga S F and Plenio M B (2007) Phys. Rev. Lett. 98 170601 

[9] Viola L, Fortunate E M, Lloyd S, Tseng C H and Cory D G (2000) Phys. Rev. Lett. 84 5466 

[10] Witthaut D, Trimborn F and Wimberger S (2008) Phys. Rev. Lett. 101 200402 

[11] Witthaut D, Trimborn F and Wimberger S (2009) Phys. Rev. A 79 033621 

[12] Witthaut D, Trimborn F, Hennig H, Kordas G, Geisel T and Wimberger S (2011) Phys. Rev. A 
83 063608 

[13] Trimborn F, Witthaut D, Hennig H, Kordas G, Geisel T and Wimberger S (2011) 

Eur. Phys. J. D 63 63 

[14] Gross M and Haroche S (1982) Phys. Rep. 93 301 (1992) Phys. Rev. Lett. 68 580 

[15] Dicke R H (1954) Phys. Rev. 93 99 

[16] Garraway B M (2011) Phil. Trans. R. Soc. A 369 1137 

[17] Baumann K, Guerlin C, Brennecke F and Esslinger T (2010) Nature (London) 464 1301 



10 



[18] Marchesoni F, Gammaitoni L and Bulsara A R (1996) Phys. Rev. Lett. 76 2609 
[19] Dalibard J, Castin Y and M0lmer K (1992) Phys. Rev. Lett. 68 580 

[20] Carmichael H J (1993) An Open Systems Approach to Quantum Optics (Springer, Berlin) 



11 



